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Objective 

• Evaluate turbulence models for integrated aircraft components 
such as the forebody, wing, inlet, diffuser, nozzle, and afterbody 

Approach 

• Integrate turbulence models into existing Navier-Stokes program 
maintaining zonal philosophy 

• Introduce corrections to baseline turbulence models to account 
for additional affects such as compressibility or separation 

• Develop algorithmic improvements for better numerical stability 
and robustness 

• Compare the strengths and weaknesses of turbulence models 

• Determine applicability of algebraic, one-equation, and two-equation 
turbulence models for typical complex flows and geometries 
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Turbulence Modeling Capabilities 

• Algebraic Models 

- Cebeci-Smith boundary layer model 

- Baldwin-Lomax boundary layer model 

- P. D. Thomas shear layer model 

• One-Equation Models 

- Baldwin-Barth 

- Spalart-Allmaras 

• Two-Equation Models 

- High Reynolds number k — e 

- Low Reynolds number k — e (Jones-Launder, Speziale, Chien, 
Lam-Bremhorst, So, and Huang-Coakley) 

- Wilcox k — <jj 

- Menter baseline and shear-stress transport blended k - u/k - e 


Navier— Stokes Time— Dependent Algorithm 

NASTD 


• Euler /Navier-Stokes Equations 

- L ami nar or Turbulent 

- Ideal Gas, Thermally Perfect Air, Equilibrium or Nonequilibrium Chemistry 

• Finite Volume Formulation 

- Roe and Coakley Flux Difference Split Schemes, Optional TVD Schemes 

• Solution Update Procedure 

- Approximate Factorization 

- Runge-Kutta Time Stepping 

- Iterative Space Marching (PNS) 

• Geometric Capabilities/Generalizations 

- Zonal Capabilities and Flexible Boundary Conditions 

- Grid Sequencing 

- Overlapping Grids 

• Turbulence Models 

- Cebeci-Smith, Baldwin-Lomax and P. D. Thomas Algebraic Models 

- Baldwin-Barth and Spalart Allmaras One-Equation Models 

- Six Low Reynolds Number k - € Models 

- k — u and Menter blended k — u)jk — t Models 
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Selected Applications 


• Transonic Supercritical Airfoil 

• Three-Element High-Lift System 

• Single Slot 2-D Ejector Nozzle 

• Confluent Mixer 

• Highly Offset 3-D Diffuser 


Modifications to Production Term 


Default calculation of production: 

p _ [l (din dujf _ 2 !duk\ 2 ' 

Re 2\dx J dxj 3 

Vorticity used in production: 


Production limiter used: 




Pj.* — min(P*., 20£> k ) = min (P k , 20 c 2 p k Re) 
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Effect of Production Limiter for the Chien k-e Model 

RAE Airfoil Analysis, Turbulent Viscosity Contours 
Mach = 0.725, a = 2.55 deg.. Re = 6.5 Million 



RAE Airfoil Analysis 

Moo = 0.725, a = 2.55°, Re = 6.5 Million 

Effect of Freestream Turbulence Level on Surface Pressure 
Cbien k — e Turbulence Model 



X 
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RAE Airfoil Analysis 

Moo - 0.725, a = 2.55°, Re = 6.5 Million 

Production Limiter Used 
Chien k — e Turbulence Model 




X 
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RAE Airfoil Analysis, Turbulent Viscosity Contours 
Mach = 0.725, a =2.55 deg., Re = 6.5 Million 
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RAE Airfoil Analysis, Mach Contours 
Mach = 0.725, a =2.55 deg.. Re = 6.5 MiUion 




NASTD Solutionof MDA Three-Element High-Lift System 


M « 0.2, AOA * 16-21 








NASTD Solution of MDA Three-Element High-Lift System 



More Accurate Solutions Have Been Obtained With One-Equation Spalart- 

Allmaras Turbulence Model 


Skin Friction Coefficients on the upper Surfaces 
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Four-Zone Grid for an Ejector Nozzle 



Primary Nozzle 


Symmetry Plane 


Single Slot Ejector Analysis 

NPR=14., Pts/Ptp=.34 

Mach Number Contours from Several Turbulence Models 
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Single Slot Ejector Analysis 

NPR=14., Pts/Ptp=.34 

Eddy Viscosity from Several Turbulence Models 






Single-Slot Ejector Nozzle Analysis 

NPR = 14 , P t jP tp = 0.34 

h!h - 100 


Comparison of Predicted Ejector Flow Rates 


Model 

wyw P 

% Error 

Experiment 

0.1010 

— 

Thomas /Baldwin-Lomax 

0.1108 

+9.7 

Baldwin-Barth 

0.1129 

+ 11.8 

Spalart-Allmaras 

0.1146 

+13.5 

Chien k — e 

0.1168 

+15.6 

Jones-Launder k — e 

0.1126 

+11.5 

Speziale k — e 

0.1127 

+11.6 

i 

So k — e 

0.1148 

+13.7 

Huang- Coakley k — e 

0.1112 

+10.1 
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Single Slot Ejector Nozzle 

Surface Static Pressure Comparison with Experimental Data 
NPR = 14.0, PJP tp = 0.34 




Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 

Surface Pressure and Computational Mesh 
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Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 
Comparison of Engine Face Total Pressures 



Pt/Pto 

I 

.95 K 



.85 

.80 


Offset Diffuser Analysis 

Ae/At=1.6, L/D=4.5, Design Pressure Ratio 
Comparison of Engine Face Total Pressures 
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Offset Diffuser Analysis 

Lower Centerline Surface Static Pressure 
A e /A t = 1.6, L/D = 4.5, Design Pressure Ratio 



0 Lower Surface Data 
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Offset Diffuser Analysis 

Upper Centerline Surface Static Pressure 
A e /A t = 1.6, LjD = 4.5, Design Pressure Ratio 


BL 

bb 

SA 

0 Upper Surface Data 







Three-Dimensional Highly Offset Diffuser 

A e /A t = 1.6, LjD = 4.5, Design Pressure Ratio 
Comparison of Engine Face Parameters 


Model 

P ’■ Uvg/Ptoo 

PtmiJPoo 

~~P _~p 

^mar ^*n:n 

P 

Experiment 

0.958 

0.890 

0.114 ; 

i 

Baldwin-Lomax 

0.936 

0.708 

0.292 

Baldwin-Barth 

0.944 

0.735 

0.265 

Spalart-Allmaras 

0.955 

0.860 

o 
i— 1 

o 

Chien k — e 

0.970 

0.894 

0.106 

Jones-Launder k — € 

0.966 

0.896 

0.104 

So k — 6 

0.975 

j 0.888 

0.112 


r. R Confluent Mixer 
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Centerline Eddy Viscosity Contours 







Comparison of Throat Total Temperatures 

BL=2.5 BL=1.5 BL=0.5 BL=0.5 BL=1.5 BL=2.5 



GE Slot Cooled Nozzle, Confluent Mixer 

Surface Temperature Distributions, 

TGX = (T/ - T IeoU )/ (T Ihot - T IaU ) 



Baldwin-Barth Mode! 

k-Epsilon Model 

Spalart-Allamaras 

CD Experiment (Average) 
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GE Slot Cooled Nozzle, Confluent Mixer 

Surface Temperature Distributions, 

TGX = (TV - T,J)/(T, U - T,^) 





Top Surface, BL= 1.5 ® e 


Bottom Surface,BL= 1 . 5 ° 
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Summary of Turbulence Modeling 
at McDonnell Dougalas Aerospace 

• The one-equation models have replaced the algebraic models as the 
baseline turbulence models. 

• The Spalart-Allmaras one-equation model consistently performs bet- 
ter than the Baldwin-Barth model, particularly in the log-layer and 
free shear layers. Also, the Spalart-Allmaras model in not grid de- 
pendent like the Baldwin-Barth model. 

• No general turbulence model exists for all engineering applications. 

• The Spalart-Allmaras one-equation model and the Chien k — e mod- 
els are the preferred turbulence models. 

• Although the two-equation models often better predict the fiowfield, 
they may take from two to five times the CPU time. 

• Future directions are in further benchmarking the Menter blended 
k — io/k — e and algorithmic improvements to reduce CPU time of 
two-equation model. 
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